Inclusive Pion Double-Charge-Exchange Above .5 

GeV 



M.J. Vicente Vacas 

Departamento de Fi'sica Teorica and IFIC, Centre mixto Universidad de Valencia-CSIC, 

E-46100, Burjassot, Valencia, Spain 

M.Kh. Khankhasayev and S.G. Mashnik 

Joint Institute for Nuclear Research, Bogoliubov Laboratory of Theoretical Physics, 
141980 Dubna, Moscow Region, Russia 

February 9, 2008 



Abstract 

A cascade model has been developed to study pion induced multichannel 
reactions ( quasielastic, SCX, DCX, absorption and vr-production) at pion 
energies above 0.5 GeV. The inclusive pion double-charge-exchange (DCX) 
reaction on ^^O,"^^ Ca and '^^^Pb nuclei in the energy range from 0.4 to 1.2 
GeV is analyzed. Pion energy spectra and double differential cross section are 
calculated. The pion production is a determinant feature in the high energy 
pion nucleus reactions, and non pion production, DCX signal is sizeable only 
at forward angles and for high energy outgoing pious. It is shown that the 
contribution to inclusive DCX processes of the conventional mechanism, with 
two (or more) quasielastic SCX steps decreases very fast as a function of the 
energy and reaches very low values at energies above .6 GeV. This opens the 
opportunity of having sizeable contributions of exotic mechanisms that are 
negligible at the A-resonance energies. 







1 Introduction 



The pion induced Double Charge Exchange (DCX) reaction has been extensively 
studied 0, 0, |[ at energies below and around the A resonance. However, many- 
questions remain not sufficiently understood, hke the angular dependence of the 
analog transitions at resonance energy 0, or the two-peaks structure of the spectra 
at forward angles that appear in the inclusive DCX process in light nuclei P, |^, 
^. There are many reasons for these difficulties. One is the existence of several 
mechanisms, some of them important, that interfere, and whose evaluation is not 
simple. Another one is the strong distortion of the pion waves at energies around 
the A resonance. Finally, for exclusive reactions, DCX processes are very sensitive 
to small details of the nuclear structure. 

At higher energies the analysis of DCX reactions is much simpler, mainly, be- 
cause the ttN cross section is much smaller, reducing the importance of distortion. 
Additionally, the contribution of some of the mechanisms, like DINT and the 
absorption |jTO| mechanism, is expected to decrease rapidly with energy. 

Furthermore, the angular dependence of the single charge exchange vrA^ reaction 
is strongly energy dependent. This has the consequence that at certain energies 
the sequential contribution to the DCX cross section results in deep minima. This 



result, found for the exclusive reaction in |11] will be shown here to be a general 
feature of inclusive DCX for all nuclei. The energies, at which the sequential mech- 
anism is very small, and distortion is probably negligible, offer the best ground for 
the investigation of some exotic mechanisms, like those involving meson exchange 



currents [jT2|, six-quark bags [0, and others involving more than two nucleons 



H^, which, although small around the A resonance, are not expected to decrease 



significantly with energy. 

We know, from our experience at low energies, how strong is the dependence 
of the DCX cross section on the nuclear structure for exclusive processes. This 
dependence is smeared out in the inclusive case, when a sum over all possible 
final states is done. Additionally, inclusive DCX offers the advantage of a much 
higher yield than exclusive experiments and does not require such a good energy 
resolution. 

At present there are no data on inclusive DCX in the energy range above .5 GeV , 
although some measurements are in progress now at ITEP (Moscow). These 
experiments have stimulated the present work. There are, however, some data 
around .5 GeV which will be used to test our model. 

The calculation presented here is a Monte Carlo simulation of the reaction. 
This is most appropriate for a multichannel situation (i.e. absorption, quasielastic, 
sex, DCX, TT-production) that renders a full quantum calculation infeasible. 

Furthermore, similar calculations [|^ || describe quite well all vr-nucleus inclusive 
reactions around the A-resonance: absorption, quasielastic scattering, single charge 
exchange, and double charge exchange. At higher energies the wavelength of the 
pion is shorter and therefore quantum interference effects should be less important. 

This paper is organized as follows: Section 2 presents the basic elements of 
the variant of the cascade model considered here (details are given in Appendix 
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A). In Section 3 the results of calculations of pion spectra (energy and angular 
distributions) for the inclusive DCX are presented, and in Section 4 the obtained 
results are discussed. 



2 Model 



2.1 Basic considerations 

A pion travelling inside a nucleus can be absorbed, can change direction, energy, 
charge, or even produce more pions. The basic inputs for our simulation will be 
the probabilities per unit length for each of these channels to happen. How these 
probabilities are obtained is presented below, in sections 2.2 - 2.4. Details on the 
simulation will be presented in Appendix A. 



2.2 TiN — > ttN 

The probability per unit length of quasielastic scattering (or single charge exchange) 
is given by 

Pn{7V^,7V^')N' = ^N{7T^,7T^')N' PN (l) 

where stands for neutron or proton, is the density of nucleons of type A^, 
and cr^(^A ^.x')]^/ is the elementary cross section for the reaction vr'*' + A^ vr'^ + A^' 
obtained from Arndt's phase shifts [ffH]. The density of protons is taken from 



experiment, and the density of neutrons is taken proportional to the density of 
protons in all results presented in this work. 

When according with the probabilities of eq. (1), a quasielastic scattering took 
place, we executed the following algorithm. First, we chose randomly a nucleon, 
of the type A^, from the fermi sea, then we boosted the vr and A^ to their center 
of mass system. Finally, we selected the scattering angle (and therefore energy) of 



the outgoing particles using again the experimental cross sections |T6[ , and boosted 
the momenta to the lab. system. When the momentum of the outgoing nucleon in 
the lab. system is below the fermi level, we consider the event to be Pauli-blocked 
and therefore keep the pion initial charge and momentum unchanged. 



2.3 Pion absorption 

Even if pion absorption is a relatively small effect at high energies (which one could 
suggest from the rapid decrease of the pion-deuteron absorption cross section at 



the energy range from 0.3 — 1.0 GeF |T7|, |T8|), there is a large number of pions 



at lower energies which are generated both by the quasielastic rescatterings and 
the pion production on the nuclear nucleons. The proportion of these pions that 
eventually comes out of the nucleus is essentially determined by the absorption 
strength. 
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Although pion absorption has been extensively studied at energies below 0.3 
GeV, little theoretical work has been done about pion absorption by complex nuclei 



at high energies, and very little experimental information is available. In ref. 
the effect of pion absorption on the pion-nucleus elastic and the SCX scattering 
has been studied in the energy range of 250 to 650 Me V". The absorptive part of 
the pion-nucleus optical potential was calculated within the framework of a many- 
body field theoretical approach. The model contains both two- and three-nucleon 
absorption mechanisms and it has been shown to agree quite well to the more 
complex microscopical model of ref. in the A resonance region. Results show a 
quite weak absorption at high energies, as was expected. Another interesting result 
presented in ref. [0] is that, whereas in the resonance region three-body- absorption 
becomes comparable with the two-nucleon mechanism, as the pion energy increases 
the effects of three-body absorption decrease again, and the two-body mechanism 
becomes dominant, as is the case at low energy. 

The probability per unit length of a pion to be absorbed is expressed in terms 
the imaginary part of the pion self-energy [| by the equation Pabs = —InT'Ilabs{k)/k. 
The imaginary part of the pion self-energy, related to two-nucleon pion absorption. 



which has been calculated in |T9[ is of the form 



ImU^lik) = -D,— p' (2) 

q 

Here, D2 = 0.0116 /m^ m6~^, k is the pion momentum in the lab system, s is 
the square of the cm. energy of a pion of momentum k and a nucleon at rest, 
a = ((73/2 + o'i/2)/3 is the spin-isospin averaged unpolarized vrA^ cross section, 
and the momentum of a virtual pion q that appears in the model of ref. |]19| is 
determined as 



^{k^ + 2mf - P 
2(F + 2m) 



which in the nonrelativistic limit actually used in [19| goes to 



q= [m(A;0-A?V2m)]i/2 

where m is the mass of a nucleon. 

The three-nucleon pion absorption has been calculated 1 19] in a similar fashion, 
and the contribution to the selfenergy is given by 



ImYi^lik) = -D.^-^ (3) 



with 



[!^(A;0-p/2m)]i/2 



^The pion selfenergy 11 is related to the optical potential as = n/2fc*' 
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and s',p' have the same meaning as s and p, but are evaluated at a kinetic energy 
of the pion equal to two thirds of the real one (T^ = |T^), and = 1.15 ■ 10""^ 
fm^mb-^MeV-^ 

The pion selfenergy pieces of eqs. (2) and (3) can readily be translated to a 
probability per unit length by the relation Pabs = —I^^abs{k)/k. 

In (see, also Erratum in p2|), the pion absorption effect on the pion- nucleus 
scattering at 800 MeV has been estimated using the quasi-deuteron (qd) model. 
Taking into account that at energies well above the resonance region the two-body 
absorption becomes dominant |]T^, it is interesting to compare the outlined above 
model (Eq.(2)) with the quasi-deuteron model. 

The imaginary part of the pion selfenergy Imllahs{k) within the framework of 
the quasi-deuteron model is given by 

/mnS(fc) =-8nTImB^2NiE)pnPp (4) 

where pp and pn are the densities of protons and neutrons. The kinematical factor 
r is determined as 

r = 7 

Pnd 

where, Ai^rA and /ijrd are the vr-nucleus and vr-deuteron reduced masses correspond- 
ingly, and 7 is the relativistic transformation factor of the tt — 2N scattering am- 
plitude from the vr-nucleus to the tt — 2N c.m.s. 

7 = [uMM^^')E2N{t^)E2N{'^')/uJ^{q)uJn{q')E2N{P)E2N{P')Y^^ 

Here, q and (f are the pion momenta before and after the collision in the vr-nucleus 
c.m.s.; K and k' are the pion momenta before and after the collision in the vr — 2N 
c.m.s., and P and P' are the total momenta of the 2N subsystem before and after 
the collision in the vr-nucleus c.m.s. We use here the "frozen" approximation which 
means that P = —2q/A and P' = P — A, where A = — g is the momentum 
transfer. 

The imaginary part of the absorption parameter B is given by 

ImB = (l/4vr)H^(T^)/(2prf(0)) 

where Pd(0) is the deuteron density at r = 0. So as in this quantity is calculated 
using the square-well potential model of a deuteron. The parameter W is related 
to the pion-deuteron total cross section as W = gcr(vr+c/ pp). The total pion- 
deuteron absorption cross section at the energy range from 0.3 GeV to 1 GeV can be 



obtained from the differential cross section data of [T^, 0. The energy dependence 
of W{T) can be approximated as 

Win) = ai/T^ + a2/Tl + a^^/T^ 

where T^^ is the pion kinetic energy in the laboratory system (measured in fm), and 
the parameters a are 

ai = 0.171, a2 = -0.612fm-\ = l.TSO/m-^ 
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2.4 Pion production 

Pion production is a determinant feature in the high energy pion nucleus reactions. 
The inelastic channels have a cross section comparable, or even larger than the 
elastic channels at energies above 0.6 GeV. Although two (or more) pion production 
channels are possible at the energies considered, the inelastic cross section is clearly 



dominated by the single pion production [23| . In this work the multipion production 
channels have been ignored 0. 

Whereas at low energies a considerable amount of data is available for most 
isospin channels, including differential cross sections, our knowledge is more frag- 
mentary in the energy regime addressed here. Data has been taken from compi- 



lation [^] and 12^ , P5[| , to obtain parametrizations of the vrA^ — ttttA^ total cross 
sections. Then, for each channel, the probability per unit length is given by the 
equation 

(5) 

When, according to this probability, a pion production event of a certain isospin 
channel has taken place, we proceed in the following way. First, a nucleon of 
the type N is randomly chosen from the fermi sea. Then the scattering angles 
and energies of the outgoing particles are selected, using the 3-body phase space 
distribution. When the momentum of the outgoing nucleon in the lab system is 
below the fermi sea level, the event is considered to be Pauli-blocked and therefore 
the pion initial charge and momentum are unchanged, and no any new pions are 
produced. 

There are five independent pion production channels induced by charged pions, 
namely. 



(1) 


tt" 


^ +p - 


TT TT n 


(2) 


TT 


+ +p - 


TT^TT^P 


(3) 




^ + n- 


^ TT^TT^n 


(4) 


TT 


+ + n - 


tt\% 


(5) 




^ + n- 


TT^TT^P 


(2) 


are 


parametrized as 



a(i) = ai„ 0.2 [1 - 0.05 (T^ - T^th)l 

cx(2) = fTi„ 0.8 [1 - 0.05 (T^ - T^th)i 

where is the inelastic tt^ + p cross section obtained from Arndt's phase shifts, 
is the kinetic energy of a pion in units of fmT^ and T^t/i is the kinetic energy 
of the pion production threshold in the same units. 



^ Note that this omission will not affect practically the higher energy part of the pion spectra, 
because for these channels some energy, at least two pion masses, has necessarily been spent. 
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In a similar fashion the other three channels have been parametrized. The 
formulas are 

o-^s) = {<^in - o-(4) - cr(5)) [1 - 0.075(r^ - T^th)], 
cx(4) = [0.30 - 5.069 x lO'X + 5.004 x IQ-^T^] 

= (^in [-2.78 X 10"^ + 0.315T^ - 4.154 x IQ-^T^] 

where is now the inelastic vr^ + n cross section, and T^r, T^^th have the same 
meaning as above. In all these cases the reproduction of data is quite satisfactory. 

Apart from the five additional channels induced by a vr" that will be obtained 
by charge symmetry, there are three channels induced by a 7r°: 



(a) 


7r° + A^- 




(b) 


7r° + iV- 


^ TT + TT^iV 


(c) 


77° + - 





where is a nucleon, tt"^ is a charged pion, and the charge of the nucleon A^' 
is determined by the charge balance in the channel (c). We use a^a) = = 
c"(c) = o'in/^, where ain is the vr^A^ inelastic cross section. These channels are less 
important, because the experiment will begin with a charged pion. Therefore the 
neutral pions are secondary pions. That means that there are fewer of them, and 
also that in the average they have less energy and will not affect much the higher 
energy part of the spectra of interest here. 

3 Results 

There are some data of pion-nucleus scattering at energies around 0.5 GeV. We 
compare the results of our program with these data, because data at higher energies 
are yet very preliminary. We should remark however that we expect our results 
to be in better agreement with experiment at higher energies, where the use of a 
semiclassical approximation like this one is more justifiable, and where the pion 
interaction with nucleons is weaker. 

In Fig.l we show the results of the present model (solid line) for the quasielastic 
7r+ scattering in ^^C, compared with the experimental data from ref. ||26|| . 

The quasielastic peak is well reproduced, in both size and width, at all angles. 
The size, being absorption of very little importance at this energy in our model, 
is governed by the elementary irN cross sections and by the fermi motion of the 
nucleons. The calculation we present has been done assuming a fermi momentum 
of 250 MeV all over the nucleus. In fact, the use of a local fermi momentum, 
obtained from the local nucleon density, results in a worse agreement with the 
data, producing narrower peaks and overestimating clearly the cross section at 
forward angles. 
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Our results underestimate the cross section at pion energies below the quasielas- 
tic peak. This seems to be a fact common to other cascade codes (see, discussion 
in ref. [^] and Sect. 4) 

Of course we cannot reproduce the elastic peak, important at low angles, given 
that we do not have collisions of the pions with the nucleus as a whole. Only pure 
quasielastic scattering is our aim here. 

In Fig. 2 we show the results of the present model (solid line) for single charge 
exchange scattering at 500 MeV. This is a very important channel for us, because 
double charge exchange requires two single charge exchange scatterings. Again, 
the quasielastic peak size and width are in a quite good agreement with data. This 
should be expected, as this peak is dominated by a single pion-nucleon collision, as 
it was the case for quasielastic scattering. The only additional information is that 
the ratio between different scattering channels is the same as in the elementary 
ttN collisions. Note also that below the peak, one can observe the same behavior 
found in quasielastic scattering. There are some pions missing in that region. In 
this figure, we show separately the pions coming from vr production, which are 
not enough to agree with the data from ||2^. Even the total suppression of pion 



absorption would not be enough to improve significantly the agreement with data. 
In Fig. 2 we also compare the data with the cascade-exiton model (GEM). This 



model of nuclear reactions |^ was proposed initially to describe nucleon induced 
reactions at bombarding energies below or at ~ 100 MeV and was later developed 
for a larger interval of bombarding energies and for the analysis of pion-nucleus 



reactions (see, e.g. and references therein.) The results of GEM, shown by 
the short-dashed curves, agree with the results of the present model. The GEM 
calculations also underestimate the low energy part of the pion spectra. 

Finally, before concentrating in the DGX channel, we will present some results 
for pion-nucleus reactions at higher energies. Our purpose is to identify the main 
features of these processes. In Figs. 3-5 we analyze the reaction vr"*" +'^^Ca —^n + X, 
and Figs. 6, 7 deal with the same reaction in lead. In all figures we split the total 
cross section into two pieces: a quasifree piece, given by those pions that come out 
of the nucleus after having only quasifree scatterings, and a pion production piece, 
given by the pions coming from events in which at least a pion production took 
place. Note that in this latter case, quasifree scatterings, prior or subsequent to 
the pion production itself, could have occurred. 

Fig. 3 shows the energy spectra of the outgoing pions. Although the energy 
spectra begin at MeV, we should remember that the results at very low energies 
are not meaningful, because the cascade method is not appropriate there. In ref. 
1^ it has been shown that quantum calculations begin to differ appreciably from 
the present kind of approach, and for inclusive processes at energies around 100 
MeV. 

Let us begin discussing the quasielastic channel. In it, we can separate a region 
of high energies, where only pions coming from one (or several) quasifree scatterings 
contribute, and a second region dominated by pions coming from vr production 
events. Note the little dip around 150 MeV in the "quasifree" part, and also the 
change of curvature in the same region of the "production" part. Both are due 



7 



to the strong absorption in the A resonance region. At energies above 300 MeV 
absorption effects are small for both the two models discussed before. 

The situation is similar for the SCX channel, although the quasifree peak is 
smaller. As expected, because DCX requires at least two scatterings, the quasifree 
peak is much smaller in its case. 

The importance of pion production channels contribution to DCX has already 
been shown at lower beam energies p[. At 600 MeV and above, DCX is totally 
dominated by vr production. Of course, this is true except for the small region 
of phase space where tt production is forbidden. Therefore, if we want to learn 
something about alternative DCX mechanisms, we should concentrate our attention 
in this region. Otherwise, any signal will be blurred by the large number of pions 
coming from production channels. And these production channels are not well 
known, neither in size, nor in the angle-energy dependence of its cross sections. 

The angular behavior of the reaction is shown in Fig. 4. All channels are 
forward peaked, essentially the "quasifree" part. Note that one of the uncertainties 
in our model is, as mentioned above, the angle-energy structure of the N{7r, 2'k)N' 
amplitudes, that has been included as being proportional to phase space. Thus, 
the "production" part reflects mostly the boost of the isotropic process from the 
CM to the lab system. 

We have selected the DCX channel to show a double differential cross section 
in Fig. 5. Given the angular behavior observed previously it is only logical to find 
that the best place to isolate a clean, non-production, DCX signal occurs at forward 
angles and for high energy pions. The rest of the spectra are totally dominated by 
the TT-production channels. Observe again the dip produced by the absorption at 
energies around the A resonance. 

Fig. 6 shows the pion spectra for the reaction vr^+^^^Pfe — vr + X at 1200 MeV. 
Many of the main features found in Calcium at 600 MeV are also relevant for this 
case. Let us mention that there is an even stronger dominance of the production 
channels, except for the quasielastic channel at the very high energy region. The 
reason is the smallness of the ttN charge exchange cross section at this energy. 
Obviously, this implies a small "quasifree" SCX cross section and an even smaller 
"quasifree" DCX cross section. 

The angular distribution, not presented here, is similar to that of Calcium, 
although a bit more forward peaked, essentially in the quasielastic channel. Fig. 7 
shows the DCX double differential cross section for the reaction tt^ +'^^^Ph -k+X 
at 1200 MeV at the same angles presented before for Calcium. Quasifree DCX is 
hardly visible with the scale used in the figure, as it could be expected from Fig. 
6. 

Finally, we have selected as observable the integrated DCX cross section, putting 
as a cut that the energy of the final pion is, at most, 150 MeV below the beam en- 
ergy. This eliminates practically all cases in which there is a pion production. The 
results, in Calcium, and as a function of the beam kinetic energy, are presented in 
Fig. 8. Very similar results are obtained putting an additional cut in angles, given 
that most pions fulfilling the previous energy condition go forward. For comparison 
we also present the quasielastic and SCX case. Note that whereas the quasielastic 
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channel changes by around a factor two (from 500 to 1300 MeV), SCX loses one 
order of magnitude and DCX decreases almost 3 orders of magnitude. 

We do not consider this as a prediction of an extremely low DCX cross section 
at high energies. We do not expect, nor claim that. First, because of a "technical" 
reason, DCX is so tiny because quasifree SCX is quite small. A change of the SCX 
vrA^ cross section by a 2% of a typical scale, say the elastic vrA^ cross section, would 
mean a factor 4 for our result. In other words, the error bars are only statistical 
and do not include the errors coming from the elementary cross sections used as 
input. A second, more important point is that the curve corresponds only to the 
ingredients of the code, namely, to two consecutive quasifree single charge exchange 
vrA^ collisions. Let us them state our result in a meaningful way: The contribution 
to inclusive DCX processes of the conventional mechanism, with two (or more) 
quasielastic SCX steps decreases very fast as a function of the energy and reaches 
very low values, compared with the quasielastic channel, at energies above 600 
MeV. 



Discussion 



In the present paper a cascade model has been developed to study pion induced 
multichannel reactions (quasielastic, SCX, DCX, absorption and vr-production) at 
pion energies above .5 GeV. 



The model has been checked by comparing with the experimental data p6|, ^ 



for quasielastic pion scattering and single charge exchange vr" scattering in C at 
Q.bGeV . We also compared the results of the present model with the CEM p8| , p9[| 
calculations. 

It has been shown that for quasielastic scattering the quasielastic peak is well 
reproduced, in both size and width, at all angles, but our results clearly underes- 
timate the cross section at pion energies below the quasielastic peak. This seems 
to be a fact common to other cascade codes as it is remarked in ref , indicating 
possibly, some interesting piece of physics missing in our description of the pion- 
nucleus reactions [^]. The CEM calculations also underestimate strongly the cross 
section at this energy range. 

For the single charge exchange scattering we also observe some missing of pions 
below the quasielastic peak. 

One could ascribe the missing cross sections to several causes that should cer- 
tainly be investigated further. For instance, at lower energies, a sizeable enhance- 
ment of the (tt, 2tt) cross sections, when comparing to quasifree calculations, had 

m 



been predicted in ref. 



P^] . The effect was related 



and was later found in ref. 
to the change of the dispersion relation of the pions in the medium. Speaking in 
simple terms, the pions are attracted by the medium. Unfortunately we do not 
know so well the pion propagation properties at higher energies, and the results 
of ref. 



cannot be easily extrapolated. In ref. 
formation of a narrow a meson, with little interaction with the medium 



26|] Zumbro et al. suggest the 
that would 



leave the nucleus prior to its decay into two pions. One should mention that the 
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obvious candidate, a weaker pion absorption, it is difficult to reconcile with the 
agreement obtained in the quasielastic peak, that for different angles is situated at 
the same energy, and also with the pion absorption data in the resonance energy 
region. 

At .6 GeV and above, DCX is almost totally dominated by tt production. To 
learn something about exotic DCX mechanisms, we should concentrate our atten- 
tion in the region of phase space where vr production is forbidden. Otherwise, any 
signal will be blurred by the large number of pions coming from production chan- 
nels. Furthermore, these production channels are not well known, neither in size, 
nor in the angle-energy dependence of its cross sections. 

The calculations for the inclusive DCX at energy 1.2 GeV show that many of 
the main features found at 600 MeV are also relevant at these higher energies. 
There is an even stronger dominance of the production channels, except for the 
quasielastic channel at the very high energy region. The reason is the smallness of 
the vriV charge exchange cross section at these energies. 

There are some interesting results in the literature concerning exclusive DCX 
processes at high energies. In particular the reactions ^^0{n^,7r~)^^Ne and also 
^^C{7r~^ ,n~y^O have been studied in ref. at energies up to 1400 MeV. Their 
resulting cross sections present two deep minima at energies around 700 and 1300 
MeV. That result does not depend on nuclear structure or the specific nuclei chosen. 
It simply reflects the energy dependence of the ttN SCX amplitude. Thus one 
expect to get a similar result for the inclusive DCX process and, possibly, gaining 
in yield and requiring a less precise energy measurement of the final pion, because 
there is no need to separate clearly a given final state of the target nucleus. To 
study this point, we have selected as observable the integrated DCX cross section 
(Fig. 8), putting as a cut that the energy of the final pion is, at most, 150 MeV 
below the beam energy. This eliminates practically all cases in which there is a 
pion production. Very similar results are obtained putting an additional cut in 
angles, given that most pions fulfilling the previous energy condition go forward. 
Whereas the quasielastic channel loses less than one order of magnitude (from 
.5GeV to l.SGeV), SCX changes by a factor 1/40, and DCX decreases 3 orders 
of magnitude. 

It should be noted that, within the present approach, the contribution to in- 
clusive DCX processes is generated by the conventional mechanism, with two (or 
more) quasielastic SCX steps. The results presented in Fig. 8 show that the sequen- 
tial mechanism decreases very fast as a function of the energy and reaches very low 
values, compared with the quasielastic channel, at energies above 600 MeV. In this 
situation, it is very important to have a possibility to compare these results with 
the experimental data and to consider other mechanisms of the DCX (like MEC) 
which might do not decrease so fast at high energies. 
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A Details on the simulation 



We generate pions, of a given momentum and charge, which travel in the z direction, 
with a random impact parameter 6, obeying |6| < i?, where R is an upper bound 
for the nuclear radius. We choose it! such that p(it!) 10~^po, with po the normal 

— * 

nuclear matter density. At the beginning, the pions are placed at the point (6, Zi^), 



with Zin = —\JR^ — l&p, and then, we proceed to move them along the z direction, 
in small steps, until either the pions get out of the nucleus or interact. 

Let us assume that P(g, r. A) is the probability of interaction per unit length, at 
the point r, of a pion of momentum q and charge A. We choose an interval 51, such 
that P{q,r, X)Sl is small compared to unity. Then we generate a random number 
X G [0, 1[. We have two possibilities: 

(a) X > P 51. In this case there is no interaction, and the pion travels a distance 
51 along the direction of its momentum q. 

(b) X < P 51. In this case there is interaction. According to its respective 
weights, we decide whether it has been absorption, quasielastic scattering, charge 
exchange, or pion production. 

When it has been quasielastic, or charge exchange, we use the procedure defined 
in section 2.2 to find the new energy, and direction of the pion, and continue to 
propagate it along its new direction, checking at every step if new interactions take 
place. 

When it has been a pion production case, we choose the channel, according 
with the respective weights given by their cross sections, and select the energy and 
direction of the two final pions using the algorithm described in section 2.4. Then 
we store the data of one of the pions and keep moving the other one. 

In the case that after moving the step 51 the pion gets out of the nucleus, and 
when the pion is absorbed, we check whether there are some other pions left inside 
the nucleus (these pions would have been produced previously at some step). If 
there are some other pions, we select one of them and begin to propagate it from 
its current position. 

When there are no pions left, we store the data -energy and angles- of all pions 
that got out of the nucleus, if any, and begin again the full procedure by generating 
a new initial pion. 

Note that between interactions pions follow a straight trajectory. Thus, even 
some classical effects hke the change in direction due to the real part of the potential 
, strong and coulombian, are neglected. 

A.l Integrated cross sections 

Let N be the total number of incident pions, and let Nc be the total number of 
events of a given channel, i.e. number of cases in which there is a 7r~ in the final 
state, then, the integrated cross section for that channel is given by 
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A. 2 Differential cross sections 

As an example, we will show the way in which angular cross sections are calculated. 
Energy distributions, or double differential cross sections arc obtained in a similar 
way. To calculate differential angular cross sections we divide the cosine of the 
polar angle in A*"^ equal angular intervals. If /i is the cosine of the polar angle with 
which a pion leaves the nucleus, we associate to it the discrete value of the angular 
bin in which it falls. Thus 

where, [.] means that we take only the integer part, and 

If we get a total number of pions in the k bin, we have 

d(Jk _ p2/ 1 
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B Figure Captions 



Fig. 1. Double differential cross section for the ^^C7(7r+, tt^)X reaction at T^=500 
MeV. Continuous line: full model.The data, taken from ref. [Q, also include the 
^"^CItt^ ,7r^)X reaction. 

Fig. 2. Double differential cross section for the ^^C(7r~, 7i^)X reaction at T^=500 
MeV. Continuous line: full model, long-dashed line: 7r-production channels, short- 
dashed line CEM model Data from ref. 

Fig. 3. Calculated spectra of the '^^Ca{7T, tt')X reactions at T7r=600 MeV. Contin- 
uous line: full model, long-dashed line: vr-production channels, short-dashed line: 
quasielastic channels. 

Fig. 4. Calculated angular distributions of the '^°Ca(7r, 7t')X reactions at T^=600 
MeV. Continuous line: full model, long-dashed line: vr-production channels, short- 
dashed line: quasielastic channels. 

Fig. 5. Calculated angle-energy distributions of the ^°Ca(7r"'", 7r~)X reaction at 
T^=600 MeV. Continuous line: full model, dashed line: quasielastic channels. 

Fig. 6. Calculated spectra of the ^^^Ph{-K, ■k')X reactions at T^=1200 MeV. 
Continuous line: full model, long-dashed line: 7r-production channels, short-dashed 
line: quasielastic channels. 

Fig. 7. Calculated angle-energy distributions of the "^^^ Ph{'K^ ,t:~)X reaction at 
T,r=1200 MeV. Continuous line: full model, dashed line: quasielastic channels. 

Fig. 8. Calculated cross section of the ^°Ca(7r, t[')X reactions as a function of the 
energy. Error bars represent statistical uncertainty. 
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